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ABSTRACT 



We present new spectroscopic identifications of 12 X-shaped radio galaxies and use the spectral data to derive starburst histories 
and masses of the nuclear supermassive black holes in these galaxies. The observations were done with the 2.1-m telescope of the 
Observatorio Astronomico Nacional at San Pedro Martir, Mexico. The new spectroscopic results extend the sample of X-shaped radio 
galaxies studied with optical spectroscopy. We show that the combined sample of the X-shaped radio galaxies has statistically higher 
black-hole masses and older episodes of star formation than a control sample of canonical double-lobed radio sources with similar 
redshifts and luminosities. The data reveal enhanced star-formation activity in the X-shaped sample on the timescales expected in 
galactic mergers. We discuss the results obtained in the framework of the merger scenario. 

Key words. Galaxies: kinematics and dynamics - Galaxies: formation - Galaxies: nuclei - Black hole physics 



1. Introduction 

A significant fraction of the most powerful FRII (Fanaroff & 
Riley type II; Fanaroff & Riley|1974| l radio galaxies exhibit two 
pairs of misaligned radio lobes associated with the same parent 
galaxy ( |Leahy & Williams|19 84 ), giving these galaxies a pecu- 
liar X-shaped appearance, which is remarkably different from 
the "canonical" double-lobed morphology of radio galaxies. 
Hence, these objects are often branded 'winged' or 'X-shaped' 
radio galaxies ( Leahy & Williams|1984 |Leahy & Parma|1992] l. 
A recent merger of two supermassive black holes (SMBHs), 
where the wings constitute the relic emission of the past radio 
jets, is one of several scenarios proposed to explain this peculiar 
radio morphology (e.g., Rottman n|2001 Merritt & Ekers|2 002 



formation mechanism for FRI and FRII X-shaped galaxies via 
backfiows (Saripalli & Subrahmanyan 2009). The required du- 
ration of these backfiows is, howe ver, too close to the typical 



radio-source lifetime of 10 yr (e.g., Rottmann 2001 Lai & Rao 



2007 1, which again argues against the backflow scenario. 



Komossa|200"6l [GeFg ely & Bi ermann|2009[ lHodges- Kluck et al.| 



2010b; M ezcua et al.|2011|i. Alternative model s employ jet-axis 
reorientation (e.g., |Dennett- Thorpe et al.|2002| l, two unresolved 
active galactic nuclei (AGN; |Lal & Rao|2007[), backflow of ma- 
terial from the main lobes into the wing s (e.g.,|Leahy & Wil liams 
1984} [Cipetti et al.||2002| |Hodges-Kluck et al.||2010a| |Landt 



et~aT]|2010[ |Hodges-Kluck & Reynolds||201 1) , or non-ballistic 
precession assuming a pre-merger state of the two black holes 
(Gong 2008 Gong et aL]|201 1) . The existence of very few X- 
shaped sources with FRI edge-darkened radio lobes (Saripalli 



|et al.|2008"t |Sar ipalli & Subrahmanyan 2009 ) was initially used 
as a strong argument against the backflow scenario. However, 
the presence in most of the FRI X-shaped galaxies of double- 
double morphologies indicative of a renewal of jet activity, to- 
gether with an alignment between the main radio axis and the 
major axis of the host galaxy, led to the proposal of a common 
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Each of the models mentioned above thus has difficulties 
in explaining the entire range of the observed properties of the 
X-shaped radio galaxies (see [Gopal-Krishna et al.||2010| for a 
review). If the observed X-shaped morphology results from a 
coalescence of a pair of SMBHs, the imprints of the preced- 
ing galactic merger should be observed in the X-shaped radio 
sources in the form of early-type host galaxies ( |Toomre|1977 1, 
statistically higher black hole (BH) masses, and enhanced star- 
formation and nuclear activity (e.g., Kauffma nn & Haehnelt| 
2000; Hopkins et al. 2008| l. It would therefore be relevant to 
investigate whether BH masses and star-formation activities dif- 
fer between X-shaped radio galaxies and canonical double-lobed 
radio sources. We initiated such a study by analyzing the opti- 
cal spectra of 29 X-shaped sources and 36 radio-loud control 
sources (cf., |Mezcua et al.||20lT| l . The results revealed that all 
the X-shaped sources analyzed are hosted by early-type galaxies, 
and that the X-shaped sample has a statistically higher mean BH 
mass and older bursts of star formation than the control sample 
of canonical radio galaxies ( [Mezcua et al.|20lTj hereafter Ml 1). 
A total of 52 X-shaped radio galaxies have been spectroscop- 
ically identified so far in several studies ( |Cheung et al.[ ]2009; 
|Landt et al.|2010[ ), but not all of these observations are suitable 
for detailed spectral analysis. To increase the statistical signif- 
icance of our studies, we identified a sample of 13 X-shaped 
radio sources with bright optical nuclei suitable for optical spec- 
troscopy observations. We present here the results of the spec- 
troscopic observations and discuss them in the context of physi- 
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cal models proposed to explain the nature of the X-shaped radio 
galaxies. 

The objects selected and the optical observations of the sam- 
ple are described in Sect. 2. Procedures applied for the analy- 
sis of the measurements are explained in Sect. 3. The results 
and discussion are presented in Sects. 4 and 5, respectively. 
Throughout the paper, we assume a A CDM cosmology with 
parameters Ho = 73 km s _1 Mpc~', Q A = 0.73, and C2 m = 0.27. 

2. Observations and data reduction 

The observed sample comprises 13 X-shaped radio galaxies re- 
trieved from a list of 100 candidates (Cheung 2007). We se- 
lected all the objects with bright (itir < 19), optical counterparts 
for which neither spectroscopic redshifts had been reported nor 
spectra were available in the Sloan Digital Sky Survey (SDSS), 
and which were observable during the period of the observa- 
tions (i.e., across a range of right ascension between 07:00:00 
and 17:00:00). 

Optical spectroscopy observations of the sample were car- 
ried out on 20-23 March 2009 with the 2. 1 -meter telescope of the 
Observatorio Astronomico Nacional at San Pedro Martir (OAN- 
SMP), Baja California, Mexico. The Boiler & Chivens spectro- 
graph was tuned to the 4000 A to 8000 A range (grating 300 
1/mm), with a spectral dispersion of 4.0 ±0.3 A/pix, correspond- 
ing to 8.0 ± 1.2 A full width at half maximum (FWHM), derived 
from the FWHM of different emission lines of the arc-lamp spec- 
trum. The instrumental resolution estimated for each object from 
the FWHM of the night-sky line at 5575 A is in the range 9-14 A. 
A 2'.'5 slit oriented along right ascension was used. To calibrate 
the spectral measurements, the spectrophotometric standard star 
Feige 34 was observed two times during every night of the ob- 
serving run. The airmass of the observations ranged between 1 .0 
and 1.3. 

For each target source, two or three exposures were taken, 
with a typical exposure duration of 1800 s (see Table 1). The 
data reduction was carried out with the irafS software fol- 
lowing standard procedures. The spectra were bias-subtracted 
and corrected with dome flat- field frames. Cosmic rays were re- 
moved interactively from all images. Arc-lamp (CuHeNeAr) ex- 
posures were used for the wavelength calibration. A spline func- 
tion was fitted to determine the dispersion function (wavelength- 
to-pixel correspondence). Sky emission lines located at known 
wavelengths were removed during the calibration in wavelength. 
The absolute flux-calibration accuracy of the spectra, provided 
by the standard star calibration, is about 15%. The relative flux- 
calibration uncertainty, estimated as the difference between the 
normalized true and estimated flux-calibration curves of the 
standard star, is < 15% for all nights. 



3. Analysis of the spectra 

The optical spectra were used to identify the most prominent 
emission/absorption lines, determine the galaxy redshifts, obtain 
mass estimates for the nuclear SMBH, and reconstruct histories 
of star-formation activity in the host galaxies. 

Since most of the observed spectra of the target objects do 
not show strong emission lines, we determined the redshift from 



Table 1 : Observation log 



Name 
(1) 


m R 

(2) 


Obs. 

(3) 


date 




Exp. time 
(4) 


J0813+4347 


16.1 


2009 


Mar 


20 


2 


x 


1800 






2009 


Mar 


22 


2 


x 


1800 


J0838+3253 


16.9 


2009 


Mar 


21 


2 


X 


1800 






2009 


Mar 


23 


2 


X 


1800 


J0924+4233 


17.8 


2009 


Mar 


23 


2 


X 


1800 


J 1008 +0030 


15.8 


2009 


Mar 


20 


2 


x 


1800 






2009 


Mar 


22 


2 


x 


1800 


J1055-0707 


17.6 


2009 


Mar 


20 


3 


X 


1800 


J 1200+6 105 


18.3 


2009 


Mar 


23 


2 


X 


1800 


J 120 1-0703 


16.4 


2009 


Mar 


22 


2 


X 


1800 






2009 


Mar 


23 


2 


X 


1800 


J1258+3227" 


17.0 


2009 


Mar 


20 


2 


X 


1800 


J1351+5559 


15.1 


2009 


Mar 


22 


2 


X 


1800 






2009 


Mar 


23 


1 


X 


1800 


J1408+0225 


18.4 


2009 


Mar 


22 


3 


X 


1800 


J 1459+2903 


16.5 


2009 


Mar 


21 


2 


X 


1800 


J1537+2648 


18.3 


2009 


Mar 


22 


2 


X 


1800 






2009 


Mar 


23 


2 


X 


1800 


J1606+0000 


15.0 


2009 


Mar 


21 


2 


X 


1800 



Column designation: (1) - object name based on J2000.0 coor- 
dinates; (2) - apparent magnitude R, taken from |Cheung ( 2007] l; 
(3) observation date; (4) exposure time, in seconds. Notes: a - 
Optical misidentification in |Cheung et al.| ([2"009). The spectrum 
of this object corresponds to a K-type star and not a X-shaped 
radio galaxy. 



the Ca II H ^3968A and K ^3934A, G band ^4302A, Mg lb 
/15175A, and Na Id /I5896A absorption lines. 

We used the stellar population synthesis code STARLIGHT 
dAsari et al.|2007HCid Fernandes et al.|2"004l[2005l[20071|Mateus| 



|et al. 2006 1 to model the observed spectra. The best fit was 
obtained by constructing a linear combination of simple stel 



lar populations (SSPs) from the stellar library of |Bruzual & 



1 IRAF is distributed by the National Optical Astronomy 
Observatories operated by the Association of Universities for Research 
in Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 



|Charlot| ([2003 ) and adding a power-law component represent- 
ing the AGN continuum emission (see Fig.[T[ e.g., Leon-Tavares 
|et al.|201 1| >. 

The STARLIGHT model for the observed spectra yields an 
estimate of stellar velocity dispersion cr t , which we corrected 
by taking into account the spectral resolution of each object, 
and from which we derived the BH mass using the empirical 
M B h - 0". relation ( |Gebhardt et al.|2000[|Tremaine et al.|2002] >, 
and the light fraction, mass fraction, age, and metallicity of the 
stellar populations used in the fit. We used these parameters to 
derive starburst histories and apply Gaussian smoothing to the 
individual starburst events in order to determine the epoch of the 
most recent starburst episode (see Fig. [2}. A quality factor Q was 
derived from the x 1 of the modeled spectra in order to quantify 
the reliability of the fit. Fits with Q > 10 can be considered reli- 
able (see Ml 1 for details). 

We estimated the rest-frame continuum flux at 5100 A from 
the SDSS photometry (Wu & Liu 2004). To assess the spectral 
classification of the X-shaped radio galaxies, we measured the 
Ca II break (Cc a n) of their absorption optical spectrum ( |Landt| 
|et al.|2002[ ) as done in Ml 1 . 

Finally, we determined the dynamic age (f a ) of the high- 
surface -brightness (active) radio lobes using 6* a /v a , where # a is 
the angular size in the FIRST image, converted to linear size as- 
suming that the sources are in the plane of the sky, and v a is the 
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4000 4500 5000 5500 6000 

Wavelength (A) 

Fig. 1: STARLIGHT fit to the spectrum of the X-shaped radio 
galaxy J1201-0703. The observed spectrum is shown in black, 
the modeled in red (for illustration purposes, displaced from the 
observed spectrum). Residuals obtained after subtraction of the 
modeled spectrum from the observed one are shown in green. 




i> 0.0100 
0.0010 
0.0001 

6 7 8 9 10 
Log Age (yr) 

Fig. 2: Age vs. mass fraction (bottom) and light fraction (top) of 
the stellar populations synthesized in the fitting of the spectrum 
of the X-shaped source J1201-0703. The squares correspond to 
stellar populations with different metallicities (in the color ver- 
sion: orange 1.0 Z and red 2.5 Z Q ). The solid curve represents 
Gaussian smoothing of the mass fraction and light fraction dis- 
tributions. The resulting peaks or bursts of star formation are 
marked with a star. 



lobe advance speed. We adopted v a ~ 0.1c, as obtained from 
observations of the aging of the synchrotron radiation spectrum 
in radio galaxies (e.g., |Scheuer| 1995) |Tingay et al.|1998] >. 

Assuming that the fueling of the low-surface-brightness 
lobes of the X-shaped sources had stopped after the high-surface 
ones were activated, the dynamic age of the passive lobes t p dur- 
ing their active stage can be estimated as 



9 P ~ t a v p 



(1) 



where 9 p is the angular size of the low-surface brightness lobes, 
converted to linear size assuming that the sources are in the plane 
of the sky, and v p is their expansion speed during the inactive 
stage. We use v„ = 0.01 c in our calculation. 



4. Results 

The new data provide the first spectroscopic identifications and 
redshifts for 8 X-shaped radio galaxies (J1008+0030, J1055- 
0707, J1200+6105, J1201-0703, J1351+5559, J1408+0225, 
J1459+2903, and J1537+2648). New spectra were obtained 
for another 4 X-shaped sources (J08 13+4347, J0838+3253, 
and J0924+4233, spectroscopi cally identified in the SDSS, 
Adelman-McCarthy et al.|2007] and for J1606+0000, identified 
by Best et al.||1999| l. The difference between our redshifts and 
the ones obtained previously are 6 Z < 0.002. The spectrum of 
these 12 sources is shown in Fig. [3] With the new identifica- 
tions, the total number of X-shaped radio galaxies spectroscop- 
ically identified is increased to 60. The last object of our sam- 
ple, J1258+3227, shows the typical stellar lines of a K-type star 
(G band, H/3, Mg I, Fe I, and Na I absorption lines; Fig. |4]l, in- 
dicating that there has been a likely misidentification (as also 
suggested by |Cheung et al.|20 09). This object is excluded from 
further analysis. Combined results from the fits to the optical 
spectra, the BH mass calculations, the epochs of the most re- 
cent starburst, and the age estimates for the radio lobes are pre- 
sented in Table 2. The columns listed are: object name based on 
J2000.0 coordinates (Col. 1), other common catalog names (Col. 
2), stellar velocity dispersion (Col. 3), BH mass derived from <t* 
(Col. 4), optical luminosity of the AGN (Col. 5) and of the host 
galaxy (Col. 6), radio luminosity (Col. 7), dynamic age of the ra- 
dio lobes (Col. 8), most recent starburst age (Col. 9), the value of 
Ca II break factor (Col. 10), spectroscopic redshift (Col. 1 1), and 
quality factor of the STARLIGHT fit (Col. 12). The total (active 
+ passive lobe) age of the radio emission is given in brackets in 
Col. 8. The objects J08 13+4347, J0838+3253, and J0924+4233 
were previously analyzed in Ml 1. We report here their new val- 
ues of stellar velocity dispersion, BH mass, most recent starburst 
age, Ca II break factor, spectroscopic redshift, and quality factor 
(0 of the STARLIGHT fit derived from the new spectra. The 
new fits have higher values of Q than in previous analysis and 
are used for the combined statistical study. 

The addition of the new spectroscopical identifications to the 
X-shaped sample of Mil provides an extended X-shaped sam- 
ple of 38 sources. This extended sample is considered from now 
on for the statistical comparison against a control sample con- 
sisting of 36 radio-loud sources with redshift z<0.3, and cov- 
ering the same range of optical and radio luminosities as the 
X-shaped sample. A control subsample of elliptical radio-loud 
sources with values of Cc a n > 0.25 and that qualify as early- 
type galaxies according to color-color diagnostics (i.e., location 
above the u-r galaxy type separator in a color-color diagram; 
Strateva et al. 2001 ) is also considered for the statistical compari- 
son (selection of both samples is described in detail in Ml 1). The 
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Fig. 3: Optical spectrum and redshift of the targeted X-shaped radio sources. *The noisy spectrum of J1200+6105 results in a higher 
uncertainty of the redshift determination for this object. 



morphological classification based on visual inspection of SDSS 
images of this subsample of elliptical control sources agrees with 
the one based on the Ca II break and color-color diagram. 



4.1. Luminosity matching 

To provide a closer statistical matching between the target and 
control samples, we used tight radio and optical luminosity 
matching between the two samples. 



We derived the optical luminosities of all of the 12 X-shaped 
radio galaxies except for two (J1055-0707 and J1201-0703), 
whose rest-frame continuum flux at 5100 A cannot be estimated 
owing to their lack of SDSS photometry. These sources are ex- 
cluded from the statistical study. Figure[5]compares the radio and 
optical luminosities in the combined target and control samples 
(with the new X-shaped objects marked). All new sources fall 
in the common range of luminosities log AL^qq^ ^ 

[43.0, 46.0] 

and log vLi.4ghz e [39.0, 44.5], defined in Ml 1 as Region 0, and 
eight of them fall within a tighter range of luminosities called 
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6000 
Wavelength [A] 



Fig. 4: Spectrum of the optical counterpart of J1258+3227, 
showing stellar lines typical of a K-type star. This object is ex- 
cluded from our analysis. 
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Fig. 5: Optical continuum luminosity versus radio luminosity at 
1 .4 GHz for the X-shaped sources of Ml 1 (filled triangles), con- 
trol sources (squares), and the new X-shaped sources analyzed 
here (empty triangles). A small subregion ("Region 1") marked 
by a rectangle inside the plot is identified to provide a tighter lu- 
minosity match between the X-shaped and control samples (see 
Mil for details). 



Reg ion 1 (with log /IZ^qq^ g [43.5, 44.25] and log v^xaghz ^ 
[40.25, 42.5]; see Mil for further details of the region defini- 
tion). 

The Kolmogorov-Smirnov (KS) test applied to the optical 
and radio luminosities of the extended X-shaped sample and the 
control sample in Region gives a probability > 5% that the two 
samples are drawn from the same parent distribution. In Region 
1, the KS-test indicates that the probability that the two sam- 
ples are the same is 87% for the optical luminosity and 19% for 
the radio luminosity. Therefore, the KS tests suggest that the X- 
shaped and control sample sources are indistinguishable in terms 
of optical and radio luminosities within Region 1 . 



1.5 



A X-shaped objects (Mil) 
□ Control sample (Ml 1) 
A X-shaped (this paper) 




2.5 



Fig. 6: Color-color diagram (g-r colors versus u-g colors) for the 
10 new X-shaped sources analyzed here (empty triangles), the 
X-shaped sources from Mil (filled triangles), and the control 
sources from Ml 1 (squares). Blackline: u-r = 2.22 galaxy type 
separator from Strateva et al. (2001 1. Sources situated above this 
line are classified as early-type galaxies. 



4.2. Host type 

The new spectra present typical features of early-type galaxies, 
showing strong metal absorption lines but no optical emission 
lines, and with values of the Ca II break of Ccaii ^ 0.3 in all the 
sources, and Ccaii ^ 0.4 in 10 of them (see Table 2, Col. 10). 
These values ensure that the galaxy is dominated by the thermal 
spectrum of the host rather than the non-thermal spectrum of an 
active nucleus or a relativistic jet. The host galaxy type can be 
inferred using a color-color diagram, obtained from the SDSS 
u, g, and r-band photometry. We apply a color-color distinction 
(Strat eva et al.|[200T} to identify the type of the host galaxy in 
the 10 new X-shaped sources considered for the statistical anal- 
ysis. The combined plot presented in Fig. [6] shows that the new 
sources also qualify as early-type galaxies. 



4.3. Black hole masses 

From the individual estimates of the BH mass, we calculated the 
mean BH mass in the extended X-shaped sample, (MBH.x-shaped), 
and derived its ratio to the mean BH mass in the control sam- 
ples r xc = <M B H,x-shaped>/<M B H,controi>- The resulting ratio is 
1.49+°^ for the X-shaped/control samples and U3+°f° for the 
X-shaped/control (ellipticals) samples. These values increase to 
1.93^Q34 and 1.50^q^, respectively, in the tighter Region 1. The 
ratios obtained are similar to the ones obtained for the previ- 
ously published sample of Ml 1, for which r xc > 1 were found in 
both Regions and 1. The median ratios are 1.03^ 25 f° r tne X- 
shaped/control samples and 1.01^2?? for the X-shaped/control 
(ellipticals) samples in Region 0, and 1.03^21? and l^+y^, re- 
spectively, in the tighter Region 1 . While the difference between 
the mean and the median values of the mass ratios is significant, 
it should be noted that application of the median is best justified 
for samples featuring a small number of strong outliers, which 
is not the case for any of our samples. Furthermore, the KS-test 
applied to the BH mass distribution of the X-shaped and control 
sample gives probabilities of 85% and 96% that the two samples 
differ significantly in Regions and 1, respectively, confirming 
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Fig. 7: Histogram of the BH mass in Region for the ex- 
tended sample of X-shaped sources (top), entire sample of con- 
trol galaxies (middle), and control ellipticals (bottom). 
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Fig. 8: Histogram of the BH mass in Region 1 for the ex- 
tended sample of X-shaped sources (top), entire sample of con- 
trol galaxies (middle), and control ellipticals (bottom). 



the difference between BH masses in the two samples expressed 
by the mean values of the mass. 

The differences between the BH masses of the two samples 
are further illustrated by the distribution of BH masses in Figs.|7j- 
[8] The histograms of the mass distributions for the X-shaped and 
control samples in Regions and 1 underline the trend of the X- 
shaped objects having statistically higher BH masses than the 
control sample. About 85% of the X-shaped sources in Region 
and 83% in Region 1 have log Mbh > 8 M Q . The percentages for 
the control sample objects are 67% and 55% in Regions and 
1, respectively. These results lend further support to the merger 
scenario in which a BH coalescence could explain the higher 
BH masses obtained in the X-shaped sample and be the possible 
origin of the X-shaped morphology. 

4.4. Starbursts 

Distributions of the ages of the most recent starburst are com- 
pared in Fig.|9]for the extended X-shaped sample and the control 
samples. Both samples present a peak of starburst activity in the 



range of 10 -10 years. The X-shaped sample notably exhibits 
a more prominent spread in the starburst ages, with 50% of the 
sources having a starburst activity older than 10 8 years. The ad- 
dition of the new galaxies to the X-shaped sample of Mil not 
only confirms this trend but also reveals a secondary peak be- 
tween 10 -10 years in the distribution of starburst ages of the 
extended X-shaped sample (see Pig . |9j . The most recent peak of 
star formation, which is present in all three samples, could be 
related to the current jet a ctivity (e.g., |Rees|1989||Labiano et aL 
2008 ; Silk & Norman|2009 1. The secondary peak, which is only 



present in the extended X-shaped sample, cannot be linked to 
the active lobes as their oldest dynamic age is ~ 10 6 years (see 
Table [2] col. 8). This strong starburst activity can only possibly 
be related to an event that occurred before the active lobes were 
formed. A delay between a merger-driven starburst and AGN 
activity is shown in recent studies (e.g., Schawinski et al.[2 010; 
|Blecha et al.||201l) . This scenario can be tested by comparing 
the starburst ages to the dynamic ages of the radio lobes for both 
the X-shaped and control samples. Histograms of the logarith- 
mic ratio of the dynamic age to the most recent starburst age are 
shown in Fig. 10 The X-shaped sources tend to have older star- 
burst ages than the dynamic ages of the radio lobes, while these 
two ages are comparable in the control sample. This trend in 
the X-shaped sample is observed for both the active lobes alone 
and the active plus passive lobes (see Fig. 10 1. The mean loga- 
rithmic ratio of the extended X-shaped sample is -1.43 ± 0.23, 
while that of the control sample is -0.14 + 0.18. According to 
the KS-test, the starburst ages of the X-shaped sample and the 
control sample differ with a probability of 98% in both Regions 
and 1. Considering the subsample of only control ellipticals, 
the KS-test gives a somewhat lower probability of 95.5% of the 
X-shaped objects being different from the control sample. 

Since the dynamic ages have been derived from the linear 
sizes of the radio lobes, we check for differences in the linear 
size distributions of the X-shaped and control samples. The KS- 
test finds that the linear size distributions differ with a probabil- 
ity of 87%. This difference is confirmed by the mean values of 
the distributions, which are (0.22 + 0.02) Mpc for the X-shaped 
sample and (0.5 + 0.2) Mpc for the control sample. The median 
values are 0.18 Mpc and 0.3 Mpc, respectively. We also study 
a possible dependence of the linear size with the BH mass and 
starburst age, but no dependences are found for any of the two 
samples. 



5. Discussion 

The origin of the X-shaped morphology in X-shaped radio 
galaxies is a matter of ongoing debate. Some recent studies in the 
X-rays ( Hodge s^ Kluck et aL|2010b] l and optical (Mezcua et al. 
201 1 1 bands and some theoretical models (Gong et al .|2011| i fa^ 



vor a recent merger of two supermassive black holes as the most 
plausible scenario. Other studies ( jHodges-Kluck et al. 2010a; 
|Landt et al.|2010| i propose a backflow from the active lobes into 
the wings to explain the peculiar morphology of X-shaped radio 
galaxies. 

In this paper, we have presented the largest statistical study 
to date of the BH mass and starburst activity in X-shaped radio 
sources. The properties of a sample of 38 X-shaped radio galax- 
ies have been compared to those of a control sample of 36 radio- 
loud active nuclei, by applying stringent luminosity and color 
criteria to provide the closest match between the two samples. 
The results suggest that a galactic merger is a more likely expla- 
nation of the origin of X-shaped sources. The merger hypothe- 
sis is further supported by the finding that all of the X-shaped 
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Fig. 9: Histogram of the ages of the most recent starburst for X- 
shaped sources (top), entire sample of control galaxies (middle), 
and control ellipticals (bottom) in Region 0. 
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Fig. 10: Logarithmic ratios of the dynamic ages of the radio 
emission to the age of the most recent starburst for Region 0. 
The top row presents the ages of active lobes in X-shaped ob- 
jects, while the second row shows the sum of the ages of the 
active and passive lobes. The two bottom rows show distribu- 
tions of the age ratios in the entire sample of control galaxies 
(third row) and the control subsample of ellipticals (fourth row). 



objects studied are hosted by early-type galaxies, which are ex- 
pected to have undergone at least one major merger event over 
the course of their evolution (Too mre[l977| l. 

Histograms of the BH mass distribution for the X-shaped and 
control samples show that X-shaped radio galaxies tend to have 
higher Mbh than the control sample. The ratio of the mean Mbh 
for the two samples confirms the trend observed, with the X- 
shaped sample having a mean BH mass that is 1.93^9"^ times 
higher than the one of the control sample in a tight common 
range of radio and optical luminosities. The higher average BH 
mass in the galaxies hosting the X-shaped radio sources may 
result from a major merger event and the consequent coalescence 
of the two central black holes, which could be the origin of the 
winged morphology observed in the radio maps. 

The occurrence of a merger event in the X-shaped radio 
galaxies should be reflected in their starburst histories. The age 



of the most recent burst of star formation and the dynamic age 
of the radio lobes was thus studied for both the X-shaped and 
the control samples. The most recent episodes of starburst ac- 
tivity were found to be statistically older in the X-shaped sam- 
ple than those found for the control sample, with 50% of the 
X-shaped sources having starburst ages older than 10 s years. In 
the X-shaped sources, the most recent starburst occurred before 
the active lobes were formed. Of particular interest is the peak 
observed in the distribution of the most recent starburst ages of 
the X-shaped sample at ages of 10 9 -10 9 5 years (1-3 Gyr). This 
peak suggests that there has been enhanced star formation due to 
a merger event, which implies a time delay of 1-3 Gyr between 
the peak of starburst activity and the end of the merger. This 
timescale agrees with the time delay of ~2 Gyr found by hydro- 
dynamical simulations of galaxy mergers ( Lotz et aLp 008 ), and 
with the timescales of 0.5-2 Gyr found observationally for indi- 
vidual objects (e.g., Ta dhunter et al.|20"05 Emonts et al.|2 006). 
On the other hand, the jet dynamic age of 10° yrs derived for the 
X-shaped sources also agrees with the time delay between the 
merger event and the onset of the radio- A GN triggered activity 
found by other studies to be 10 6 years (e.g., |Emonts et al. 2006). 
These results assume a jet advance speed of 0.1c. Estimates 
of the lobe advance speed based on synchrotron aging studies 
can reach 0.2c for the most powerful radio galaxies, while the 
study of three samples of very powerful radio sources yielded 
estimates below 0.1c (with a poor statistical significance; i.e., 
Scheuer 1995). The use of an advance speed as low as 0.01c in 
our calculations would increase the dynamic age of the active 
lobes by an order of magnitude. However, this would not affect 
the results derived from the ratio of dynamic to starburst ages, 
which would still indicate that the most recent starburst activ- 
ity occurred before the active lobe formation. This result is also 
essentially unchanged when an advance speed of 0.2c is used. 

All the results obtained in our study strengthen the hypothe- 
sis that X-shaped radio galaxies were formed by mergers. More 
detailed studies of star formation and gas kinematics in these 
objects will yield stronger constraints of this scenario. 
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